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This article reports on recent searches for single-top-quark production by the CDF collabora- 
tion at the Tevatron using a data set that corresponds to an integrated luminosity of 955 pb -1 . 
Three different analyses techniques are employed, one using likelihood discriminants, one neu- 
ral networks and one matrix elements. The sensitivity to single-top production at the rate 
predicted by the standard model ranges from 2.1 to 2.6 a. While the first two analyses observe 
a deficit of single-top like events compared to the expectation, the matrix element method 
observes an excess corresponding to a background fluctuation of 2.3 a. The null results of the 
likelihood and neural network analyses translate in upper limits on the cross section of 2.6 pb 
for the i-channel production mode and 3.7 pb for the s-channel mode at the 95% C.L. The 
matrix element result corresponds to a measurement of 2.7^'g pb for the combined t- and 
s-channel single-top cross section. 

In addition, CDF has searched for non-standard model production of single-top-quarks via 
the s-channel exchange of a heavy W' boson. No signal of this process is found resulting in 
lower mass limits of 760 GeV/c 2 in case the mass of the right-handed neutrino is smaller than 
the mass of the right-handed W' or 790 GeV/c 2 in the opposite case. 
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1 Introduction 

In pp collisions at the Tevatron top quarks are mainly produced in pairs via the strong force. However, the 
standard model also predicts the production of single top-quarks by the weak interaction via the s- or t- 
channel exchange of an off-shell W boson. While early Run II searches by the CDF and D0 collaborations, 
based on data sets corresponding t o 162 pb" 1 , 230 pb" 1 or 695 pb" 1 of integrated luminosity, did not find 
evidence for single-top productionDEES one f the latest D0 analyses^ using data with L mt = lfb" 1 
observes an excess of single-top-like events of 3.4 a. 

The single-top production cross section is predicted to a s+t = 2.9 ± 0.4 pb for a top mass of 175 
GeV/c 2 ^l which is about 40% of the top-antitop pair production cross section. The main obstacle in 
finding single-top is however not the production rate of the signal but the large background rate. After 
all selection requirements are imposed, the signal to background ratio is approximately 1/16. This 
challenging, background-dominated dataset is the main motivation for using multivariate techniques. 

2 Standard Model Searches 

In this article we present three new CDF searches for standard model single-top production and one new 
search for a W boson decaying into tb. All four analyses are based on the same event selection and use 
the same Run II data set corresponding to an integrated luminosity of 955 pb -1 . The event selection 
exploits the kinematic features of the signal final state, which contains a top quark, a bottom quark, 



Table 1: Expected number of signal and background events and total number of events observed in 955 pb in 

the CDF single-top dataset. 



Process 


N events 


Process 


N events 


W + bb 
W + cc 
Wc 

Mistags 


170.9 ±50.7 

63.5 ± 19.9 

68.6 ± 19.0 
136.1 ± 19.7 


non-W 
tt 

Diboson 
Z + jets 


26.2 ± 15.9 
58.4 ± 13.5 
13.7 ± 1.9 
11.9 ±4.4 


Total background 


549.3 ± 95.2 


i-channel 


22.4 ±3.6 


s-channel 


15.4 ± 2.2 


Total prediction 


587.1 ±96.6 


Observed 


644 



and possibly additional light quark jets. To reduce multijet backgrounds, the W originating from the 
top quark is required to have decayed leptonically. One therefore demands a single, isolated high-energy 
electron or muon (Er(e) > 20 GeV, or Pt(a*) > 20 GeV/c) and large missing transverse energy from the 
undetected neutrino _^t>25 GeV. Electrons are measured in the central and in the forward calorimeter, 
1 77 1 < 2.0. To further suppress events in which no real W boson is produced, called non-W background, 
additional cuts are applied. The cuts are based on the assumption that these events do not produce $t 
by nature but due to lost or mismeasured jets. Therefore, one would expect small $t and small values 
of the angle between IpT and a jet. We further reject dilepton events from Z decays by requiring the 
dilepton mass to be outside the range: 76 GeV/c 2 < M« < 106 GeV/c 2 . 

The remaining backgrounds belong to the following categories: Wbb, Wcc, Wc, mistags (light quarks 
misidentified as heavy flavor jets), non-VF and diboson WW, WZ, and ZZ. We remove a large fraction of 
the backgrounds by demanding exactly two jets with Et > 15 GeV and \r]\ < 2.8 be present in the event. 
At least one of these two jets should be tagged as a b quark jet by using displaced vertex information from 
the silicon vertex detector (SVX) . The numbers of expected and observed events are listed in table [TJ 

2.1 Likelihood Discriminant Analysis 

One multi-variate analysis uses likelihood discriminants to combine several variables to a discriminant 
to separate single-top events from background events. One likelihood discriminant is defined for the 
f -channel, one for the s-channel search. Seven or six variables are used, respectively. The likelihood 
functions are constructed by first forming histograms of each variable. The histograms are produced 
separately for signal and several background processes. The histograms are normalized such that the 
sum of their bin contents equals 1. For one variable the different processes are combined by computing 
the ratio of signal and the sum of the background histograms. These ratios are multiplicatively combined 
to form the likelihood discriminant. 

One of the variables used in the analysis is the output of a neural net b tagger. In figure QJl the 
distribution of this b tag variable is shown for the 644 data events. In case of double-tagged events the 
leading b jet (highest in Et) is included in this distribution. The neural net b tagger gives an additional 
handle to reduce the large background components where no real b quarks are contained, mistags and 
charm-backgrounds. Both of them amount to about 50% in the W + 2 jets data sample even after 
imposing the requirement that one jet is identified by the secondary vertex tagger of CDF". 

The ^-channel likelihood function is shown in figure [TJd. The best sensitivity (expected p- value 2.5%) 
is reached by combining the two likelihood discriminants in a two-dimensional fit where ^-channel and s- 
channel are considered as one single-top signal (combined search). The observed data show no indication 
of a single-top signal and are compatible with a background-only hypothesis (p-value 58.5%). The upper 
limit on the combined single-top cross section is found to be 2.7 pb at the 95% C.L., while the expected 
limit is 2.9 pb. The best fit for the cross sections yields a t = 0.2+^2 pb and a s = O.lj^Jpb. 

2.2 Neural Network Search 

In the second analysis a neural network is used to combine 23 kinematic or event shape variables are com- 
bined to a powerful discriminant. Figure^ shows the observed data for the combined search compared 




Figure 1: (a) Output distribution of the neural net b tagger for 644 candidate events in the W + 2 jets bin. 
Overlayed are the fitted components of beauty-like, charm-like and mistag templates The yellow error band 
indicates the statistical uncertainties of the fitted sum. (b) The distributions of the t-channel likelihood function 
for CDF data compared to the Monte Carlo distribution normalized to the expected contributions. A logarithmic 

scale is used. 



to the expectation in the signal region defined by a neural network output between 0.3 and 1. In the 
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Figure 2: Single-top search with neural networks at CDF: (a) Data compared to the standard model expectation 
in the signal region (neural network outputs larger than 0.3). (b) Likelihood contours of the separate neural 

network search. 



combined search where the ratio of i-channel and s-channel cross sections is fixed to the standard model 
value a p-value of 54.6% is observed, providing no evidence for single-top production. The corresponding 
upper limit on the cross section is 2.6 pb at the 95% C.L.. To separate t- and s-channel production two 
additional networks are trained and a simulantcous fit to both discriminants is performed. The best fit 
values are at — O^tJ'Jjpb for the i-channel and a s = 0.71q ypb. The corresponding upper limits are 2.6 
pb and 3.7 pb, respectively. The observed p-value is 21.9%. 



2.3 Matrix Element Analysis 

Another way to discriminate signal from background is to compute leading order matrix elements for 
signal and background processes. The measured four-vectors of the jets and the charged lepton are 
used as experimental input. Constraints on energy and momentum conservation are applied. The jet 
energy measurements are corrected to parton level energies using transfer functions. To obtain a relative 
weight for a certain hypothesis one integrates over jet energies and the momenta of the incoming partons 
using the parton density functions as integration kernel. The weights for the individual hypotheses are 
combined event-by-event by forming a ratio signal and signal-plus-background weights. The resulting 
discriminant is named event probability density (EPD). The measured EPD distribution is shown in 
figure [31 compared to the fitted event rates for the different processes. A single-top signal corresponding 
to a 2.3 a excess is observed. The associated single-top cross section is 2.7lj;|pb. 




Figure 3: (a) Event probability density (EPD) based on matrix elements. The inset shows the signal region with 
EPD > 0.7. (b) W' — ► tb search: The invariant mass of the charged lepton, the reconstructed neutrino and the 

two leading jets. 



3 Search for a W boson 

Based on a the same event selection as the standard model searches CDF has also searched for a W 
boson in the decay channel, W — > tb. The signal is modeled using the event generator pythia^. The 
invariant mass of the charged lepton, the reconstructed neutrino and the two leading jets is used as a 
discriminant. The measured data are shown in comparison to the standard model prediction in figure^. 
No evidence for a resonant W' boson production is found, yielding limits on the production cross section 
ranging from 2.3 pb at M(W) = 300 GeV to 0.4 pb at M(W) = 950 GeV. Utilizing theoretical cross 
section calculations^, lower limits on the W' mass are set: M(W') > 760 GeV if the mass of potential 
right-handed neutrinos is below M(W) and M(W') > 790 GeV otherwise. 

4 Conclusions 

The collaboration has performed searches for standard and non-standard model single-top production 
using a data set corresponding to an integrated luminosity of 955 pb -1 . No evidence for these processes 
could be established. Two standard model searches based on likelihood discriminants or neural networks 
find no excess which can be attributed to single-top production, while the matrix element analysis finds 
an excess of 2.3 a compatible with single-top production. The overall consistency of all there analysis is 
only 1%, since the analyses feature a correlation between 60% and 70%. At present, there are no hints to 
other causes than statistical fluctuations. Even larger data samples will be available in the near future 
and will help clarify the situation. The W search in the single-top channel establishes new upper limits 
of M(W) > 760 GeV or M(W) > 790 GeV depending on the mass of the right-handed neutrino. 
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